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The quantum adder is an essential attribute of a quantum computer, just as clas- 
sical adder is needed for operation of a digital computer. We model the quantum full 
adder as a realistic complex algorithm on a large number of qubits in an Ising-spin 
quantum computer. Our results are an important step toward effective modeling 
of the quantum modular adder which is needed for Shor's and other quantum algo- 
rithms. Our full adder has the following features: (i) The near-resonant transitions 
with small detunings are completely suppressed, which allows us to decrease errors 
by several orders of magnitude and to model a 1000-qubit full adder. (We add a 
1000-bit number using 2001 spins.) (ii) We construct the full adder gates directly 
as sequences of radio- frequency pulses, rather than breaking them down into gen- 
eralized logical gates, such as Control-Not and one qubit gates. This substantially 
reduces the number of pulses needed to implement the full adder. [The maximum 
number of pulses required to add one bit (F-gate) is 15]. (iii) Full adder is realized in 
a homogeneous spin chain, (iv) The phase error is minimized: the F-gates generate 
approximately the same phase for different states of the superposition, (v) Modeling 
of the full adder is performed using quantum maps instead of differential equations. 
This allows us to reduce the calculation time to a reasonable value. 

PACS numbers: 03.67.Lx, 75.10.Jm 
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I. INTRODUCTION 

A quantum computer (QC) could efficiently solve some important problems using the 
superposition principle of quantum mechanics if the number of qubits in the register of the 
QC is sufficiently large. However, with current technologies, it is very difficult to imple- 
ment a quantum computer with many qubits. It is therefore of importance to simulate and 
test quantum algorithms on digital computers. One of the main obstacles for experimental 
implementation of the quantum information processing is decoherence caused by interac- 
tion of the quantum system with the environment. The second obstacle is inaccuracy in 
implementation of quantum protocol. In this paper we neglect these two causes of errors. 

Since an Ising spin QC based on nuclear or electron spins is operated by radio-frequency 
(r/) pulses, the wavelength of each pulse is much larger than the distance between qubits 
and much larger than the size of the whole quantum register. That is why an Ising spin 
QC is characterized by the third cause of error — nonlocality of interaction of the electro- 
magnetic waves with the qubits of the register. Since each r/ pulse affects all spins in the 
chain, the state of each qubit of the register has a small probabihty of being changed by this 
pulse. For typical Ising spin QC parameters the probabilities of the unwanted excitations 
of the qubits are very small. This problem is not very important in a conventional NMR 
spectroscopy or in a QC with a small number of qubits because the number of pulses is 
relatively small. The operation of a QC with large number of qubits requires large number 
of pulses. Consequently errors can accumulate and a theoretical approach and numerical 
simulations are required to estimate this type of error. 

The Hilbert space of a many-qubit QC is exponentially large. This feature creates two 
major problems, (i) Even if the number of initial states is small, the number of the states 
created by the pulses of a protocol increases exponentially, (ii) Due to the nonlocality of 
interaction of the r/ pulses with the qubits, the direct simulation of the dynamics requires 
either the solution of large system of coupled differential equations for a long period of time 
or the diagonalization of large matrices. In order to overcome these difficulties a perturbation 
theory was developed in our previous papers [1-9] . In this paper we apply our perturbation 
approach to simulate the dynamics of a full adder with 1000 qubits. (We suppose that 
there are / = 1000 addend qubits. The total number of qubits including carry-over qubits is 
L = 2/ + 1 = 2001.) 
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In order to solve problem (i), we formally divide all states of the quantum register in 
"useful" states and "unwanted" states. The useful states are the states which realize the 
quantum algorithm. The unwanted states are the states which are created from the useful 
states by the (unwanted) action of pulses of the protocol. For typical QC parameters the 
probabilities of generation of the unwanted states from the useful states are small (of the or- 
der of fi, where /i ~ 10~^). The probabilities of generation of the unwanted states from other 
unwanted states are of the order of /i^, so that the latter states can be neglected. In spite 
of the fact that the total number of populated unwanted states increases exponentially, the 
number of unwanted states with sufficiently large probabilities (of the order of fi) increases 
only linearly with the number of pulses. 

For solution of problem (ii), we formulate the quantum dynamics in terms of quantum 
maps. We use 

amplitudes of unwanted transitions. In this approach, one pulse of the protocol corresponds 
to one discrete step of the map. 

The full adder in the Ising spin QC was simulated in Ref. 0. Advanced features of 
the quantum full adder presented in this paper are enumerated in Abstract. In Sec. II we 
consider the dynamics of the Ising spin QC. The implementation of full adder in the Ising 
spin QC is considered in Sec. III. Quantum protocols for the full adder are described in Sec. 
IV. The quantum map approach is analyzed in Sec. V. Our numerical results are presented 
in Sec. VI. In Sec. VII we draw our conclusions. 

II. ISING SPIN QUANTUM COMPUTER 

The Hamiltonian for the Ising spin chain placed in an external permanent magnetic field 
and driven by rectangular r/ pulses can be written in the form: 

L-1 L-2 Q{n) L-1 

Hn = -Y: Y: m+i-^ E {4^exp i-t {uj + v^n)] + h.c] = H, + V^{t). (1) 

fe=0 fe=0 fc=0 

Here h = 1] = ±. /|; Jjf, I^, and are the components of the operator of the fcth 
spin 1/2; uok is the Larmor frequency of the fcth spin; J is the Ising interaction constant; 
n*-"-* is the Rabi frequency (frequency of precession around the resonant transversal field 
in the rotating frame); Un is the frequency of the pulse; and ipn is the phase of the nth 
pulse. The Hamiltonian (Q) is written for the nth rectangular r/ pulse. We assume that the 
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Larmor frequency of kth spin is Uk = Wq + kSu, so that the Larmor frequency difference 
6uj = cijfc+i— cijfc between the neighboring spins is independent of the spin number, k. Below we 
omit the index n which indicates the pulse number. The long-range dipole-dipole interaction 
is suppressed by choosing the angle between the chain and the external permanent magnetic 
field to be equal to the magic angle |1Q] . Note, that the Hamiltonian (0) allows the transitions 
associated with flip of only a single spin. 

Because of the magnetic field gradient each spin has its unique Larmor frequency. The 
selective excitations of the spins in the chain are performed by choosing the frequency of the 
r/ pulse to be equal (or approximately equal) to the Larmor frequency of the spin which we 
want to excite. 

A. Suppression of the near-resonant transitions 

The constant Ising interaction between the spins is characterized by the Ising constant 
J. Due to this interaction, the transition frequency of kth spin depends on the states of 
(A;— l)th and (A;+l)th spins, so that the same spin can have different transition frequencies for 
different quantum states of a superposition. The kth spin has different transition frequencies 
for the four possible kinds of states: 

I • • • Ofc_infclfc+i . . .), I . . . Ifc-i^fcOfc+i . . .), I . . . Ofc_infcOfc+i . . .), | . . . Ik-iUklk+i ■ ■ •), 

(2) 

where can be equal to or 1. [The first and second states in Eq. Q have the same 
transition frequencies.] The interaction between the spins is necessary for implementing the 
conditional quantum logic in the QC. 

In order to ensure optimal selective excitations, a large magnetic field gradient (of the 
order of 10^ T/m [11-13]) is required, so that the inequality J <^ 5uj is satisfied. The 
differences between the transition frequencies of kth spin in the states in Eq. Q are of 
order of J . If one wants to flip the kth spin in one of the states in Eq. Q and does 
not want to flip the same spin in the other two or three states, then one must to suppress 
these unwanted transitions (near-resonant transitions) with the detunings (from resonance 
condition) of the order of J . Since J ^ 5uj, the detunings for the near-resonant transitions 
are much smaller than the detunings for the nonresonant transitions characterized by 5uj. 
That is why the near-resonant transitions in the general case generate the largest errors. 
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The approach which allows one to completely suppress the near-resonant transitions was 
developed in Ref. Q]. In order to flip the kth qubit in the flrst two states in Eq. Q and to 
suppress the transitions in the third and fourth states, the value of the Rabi frequency fl 
should satisfy the 27rii"-condition [l^ 

^{K)= K = l,2,.... (3) 



The duration of this pulse is r = tt/Q{K) (vr-pulse) and the frequency is equal to the Larmor 
frequency = Wq + k5uj of the fcth spin. Here the upper indices of indicate the states 
of the {k — l)th and {k + l)th spins and = z/^°. The probability of the spin flip is 
independent of the phase of this pulse. We denote the pulse with these parameters as 
Q^^{ip), where the two upper indices indicate the states of the neighbors of the kth spin, 
Qk^if) = Qk^i'^)y k is the number of the spin to be flipped, and ip is the phase of the pulse. 

In order to flip the kth qubit in the third or fourth state in Eq. (j21) and to suppress the 
transitions in the other three states, two pulses are required The flrst pulse flips the 
kth spin and the second pulse removes the unwanted states created by the near-resonant 
transitions (generated by the flrst pulse) from the register of the QC. 

The parameters of the flrst pulse required to flip the kth spin in the state 
I . . . Ok-iUkOk+i ■ ■ •) have the following values the Rabi frequency is 

2 / 

Q,{K2) = ^^—, ^2 = 1,2,...; (4) 

- 1/4 



the frequency is = wq + k6uj + 2 J; the time-duration of the pulse is r = 7r/r22(-f^2); and 
the phase ip of the pulse is arbitrary. 

The parameters of the second pulse required to remove the unwanted states created by 
the flrst pulse have the following values. The Rabi frequency is 

2 / 

^c{K,)= , ^ , K, = l,2,...; (5) 



nK, 

TT+I3J 



the frequency is equal to the Larmor frequency, WQ + k6uj, of the kth spin. The time-duration, 
Tc, and the phase, ^9°°, of the pulse are: 

rc = ^^, ^T = + ^ + 2Jto + Q, (6) 

where to and (p are, respectively, the beginning time and the phase of the flrst pulse, the 
upper index of indicates that the flrst pulse is used to flip the kth spin in the state 
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I . . . Ok-iTikOk+i . • •). The new parameters introduced in Eqs. (jSJ and © depend only on the 
value of K2: 



tan + , e = tvJkI - 1/4, tan^ = , ^ sin 9. 



^i^l + 3/4 V2V 2 /y V 2 /, M ^^2_i/4 

(7) 

The parameters of the two pulses required to flip the A;th spin in the state 
I . . . Ik-iUklk+i ■ ■ ■) have the same values as those for the state | . . . Ofc-in^Ofc+i . . .), dis- 
cussed above, except that the frequency of the first pulse, u^^, and the initial phase, ipl^, of 
the second pulse are 



" - Wo + k6uj - 2J, ipl^ = -e + ip- 2 Jto - ©• (8) 



Below we treat the two pulses required to flip the kth spin in the state | . . . Ok-iUkOk+i ■ ■ ■) 
as one combined pulse and denote it as Qf{ip). The notation Ql^^^ip) is used for the two 
pulses required to flip the fcth spin in the state | . . . Ik-ink^k+i ■ ■ ■)■ 

The pulses acting on the 0th qubit and (L — l)th qubit, Ql{y:>) (here k = 0, L — 1 and 
the upper index indicates the state of the neighbor) have the following parameters: 



k 



Wo + k6u + J, T = TT/n{K), (9) 



where Q{K) is defined by Eq. Q and the phase does not affect the probability errors 
generated by this pulse. The pulse QQ^f) flips the 0th qubit in the state | . . . Oino) (^0 = 0, 1) 
and does not flip the 0th qubit in the state | . . . li^o). The pulse Q^^^ii^p) flips the {L — l)th 
qubit in the state |TiL-iOL-2) and does not flip the (L — l)th qubit in the state \nL-ilL-2)- 
The pulse Ql{ip), k = 0, L — 1, has the frequency z/^ = wo + k6uj — J and the other parameters 
are the same as for the pulse Ql{ip). 



III. QUANTUM FULL ADDER 

The quantum full adder FA, first suggested in [14], adds a number A to a superposition 
of numbers Bi which are coded by quantum states in the register of a QC, 

M M 

FA{A) Y^CsAom) = E^G.(T)|G,)e-^^«.^, (10) 

i=l 1=1 

where the state \Gi) is obtained as the result of summation of the state \Bi) and the number 
A, the addition is performed in the interaction representation jl^, Eq. = {Gi\HQ\Gi), M 



7 



a 


b 


c 


s 


c 























1 


1 








1 





1 





1 








1 








1 


1 





1 


1 





1 





1 


1 


1 








1 


1 


1 


1 


1 


1 



TABLE I: Table for binary addition. 

is the number of the states in the superposition, T is the total time-duration of the full 
adder protocol, (7^,(0) and Cg^(T) are the complex coefficients satisfying the normalization 
condition I]^=o^ IC'nl^)!^ = 1- For the ideal full adder these coefficients before and after 
implementation of the full adder are equal to each other, Cq^^^{T) = Cb,(0). The values of 
the numbers Bi and Gi are defined by the states of the spins in the spin chain. In binary 
notation, the orientation of a spin along the direction of the permanent magnetic field 
corresponds to the bit and the orientation of a spin in the opposite direction corresponds 
to the bit 1. Since the quantum logic is implemented in the quantum computer we use 
the term "qubit" (quantum bit) instead of the (classical) term "bit" for the numbers coded 
through the states of the spins. To implement the addition, one should also include the 
carry-over qubits. The addition pO|) is realized by fiipping definite spins in definite states 
[see Eqs. (dH) - (O below] in the sum in Eq. (fTIUl . 

The quantum full adder adds numbers using the same rules as a classical full adder. A 
classical full adder operates with an input of two addend bits, a and b, and a carry-over bit 
, c, as shown in Table I, where C and s are, respectively, the output carry-over and sum. 
The latter can be expressed as s = a © 6 © c, where © is addition modulo 2. The output 
carry-over is expressed as C = a6 © ac © be. 

The quantum full adder (fTUj) consists of a series of gates F{a) (F-gates), where a is the 
addend bit of the number A. There are two types of gates: F{0) adds the bit of the number 
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A, and F{1) adds the bit 1 of the number A. In other words, the value of the number A 
determines which protocol is applied to the superposition of the numbers Bi. The number 
of the F-gates is equal to the number of the addend qubits. Initially a number Bi in the 
register is represented in the form 

\B^) = |6t-'lOL-2&L-'30L-4 • • • blOsb%Oo), (11) 

where we omit the index i (which indicates the state number) in the right-hand side, the 
lower indices in the right-hand side indicate the spin number and the upper indices indicate 
the addend qubit number of the number Bi. The action of the first F-gate is defined as 

F2,i,oi(^')\...b%b%Oo) = \...b%CXs'o), (12) 

where the lower indices of F2^ifl{a^) show that this gate acts on the 0th, 1th, and 2nd spins. 
The output carry-over qubit of the first gate, 0°, is used as an input carry-over of the next 
(second) F-gate, -F4^3^2(fl^)- Using the notation introduced in Table I, this can be written in 
the form = The F-gate acts on intermediate qubits as 

Fk+i,k,k-i{a^)\ ■ ■ ■ bi^,0,4-ibi-2 . . .) = I . . . Ci^.i6i4-iC2 • • (13) 
where k = 2j + 1. 

If the addend numbers consist of I qubits, the FA (A) gate can be written in the form 

FA(A) = Fi_i,i_2,L-3(a'-') ■ • ■ i^4,3,2(«')i^2,i,o(a°), (14) 

where |A) = |a'~^a'~^ . . . a^a°), L = 2/ + 1 is the total number of spins and the gate sequence 
should be read from right to left. 

IV. QUANTUM PROTOCOLS FOR F-GATES 

The three-qubit F-gates can be implemented using the technique developed in Ref. [l| for 
one and two-qubit gates. The pulses Ql'^^ip), Ql^{ip), and Ql^{ip) (Q-pulses), considered in 
Sec. Ill A| flip the kth qubit only in the states with the definite orientations of {k — l)th and 
{k + l)th qubits. [We also mention here that Ql^{{p) = Ql^{(p).] We say that these pulses 
are probability-corrected. Using the Q-pulses for implementation of the F-gates make the 
probability-corrected F-gates. However, all these pulses in general generate different phases 
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TABLE II: Phases generated by the Q-pulses acting on intermediate qubits. The asterisk indicates 
that the resonant transition from the state shown in the first column of the table to the other state, 
associated with the flip of the A;th qubit, takes place. The phases 6, a, 9 and 7, are defined in 
Eqs. and (fT3|) . 

for different states of the superposition, and the differences between these phases for different 
states are not small. On the other hand the F-gates should generate the same phase for all 
states. In this case we say that such a protocol is phase-corrected. To make the protocols 
for F-gates phase-corrected one should choose the proper set of the phases if for the different 
Q-pulses of an F-gate. 
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TABLE III: Phases generated by the Q-pulses acting on the right edge (0th) qubit. 

The phases generated by the Q-pulses acting on different states were calculated in Ref. [l| 
for the case K2 = K and = K, when ^2 = 2Q and Qc ~ 21]. In order to decrease the 
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TABLE IV: Phases generated by the Q-pulses acting on the left edge {L — 1th) qubit. 

errors we chose smaller values for the Rabi frequencies ^2 ~ ^ ai^d ~ ^ and take 
K2 = Kc = 2K. The phases generated by these Q-pulses for different states are shown in 
Table II, where the phases 9 and O are defined in Eq. (jZj) and jl| 



a 



Tc^K^ - 1/4, 7 = ^{nK.y-i-K + py. 



(15) 



The phases generated by the Q-pulses acting on the edge qubits are shown in Tables III and 
IV. 

The phase- and probability- corrected F-gates for intermediate qubits are defined in Table 
V. The F-gates for the left addend qubit are defined in Table VI. We must note that the 
sets of phases if in the Tables V and VI are not unique since the number of linear equations 
(which is equal to^ 8) for finding the phases is smaller than the number of variables, N, in 

[A^ = 13 for F(0) and N = 15 for F(l).] The F-gates for the right 



these equations |l|. 
addend qubit are 



i^2,i,o(o) = gj(o)g^^(o)g?M-^ + 3a 



(16) 
(17) 



The gates (fT^ and (fT7|) produce the overall phase — Sa. 

We should note that our full adder protocols are formulated in terms of rf pulses but 
not in terms of generalized logical gates, such as Control-Not and one-qubit gates. Our 
approach requires much fewer pulses which allows us to decrease the errors significantly. 
For comparison, implementation of the F(0) gate requires one Control- Not gate and one 
Control-Control-Not gate. The gate F(l) also needs one more Control-Not gate and Not 



gate 



m 



. Implementation of only one Control-Not gate requires 15 or less pulses 



1. In 



11 



pulse 


-^fc+i,fc,fc-i(0) 


-Pfc+l,fc,fc-l(l) 


number 


pulse 


phase if 


pulse 


phase if 


1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 


Qf 

r)00 

Qf-i 

Qf-1 

Qf-, 

Qf 

Qf 

Qf 

Qf^ 

Qf 
Qf 


-7-26 
77 - 86 + + 14a 
1-47 + 46-30 
-6 + ha 

vr + 57 - 76 + + 5a 
TT + 27 - 26 + + 2a 


-87 + 106 -20 -11a 
TT - 57 + 76 - i0 - 7a 

e-a 






Qf 

Qf-. 

Qf-. 

Qf-. 

Qf 

Qf 

Qf 

Qf-. 

Qf 

Qf 

Qf-. 

Qf 

Qf 


IT -2-^-^^-2,0 + a 
-f + 37-6 + |0 + a 

-f -7 + 36 + i0 
-f -37 + 36- i0 + a 
-47 - 30 - a 
-f-0 




vr - 37 + 36 - i0 - a 








TABLE V: Protocols for F-gates performing addition of an intermediate addend qubit. 

this paper we use the same number of pulses for implementation of the complete F-gate 
protocols. 

V. QUANTUM MAPS 

As discussed in Introduction, an exact numerical simulation of the dynamics of the 1000- 
qubit quantum computer is practically impossible. However, as shown in this Section, simu- 
lation of the dynamics using a perturbation approach is possible. The only source of errors 
in our model of quantum computer are the nonresonant transitions. Wc use analytical for- 
mulas for the amplitudes of the nonresonant transitions to formulate the dynamics in terms 
of quantum maps. One step of the map corresponds to a single Q-pulse of the protocol. Each 
map consists of the following four successive steps: (a) Each useful state of a superposition 



12 



pulse 


-^fc+i,fc,fc-i(0) 


-Pfc+i,fc,fc-i(l) 


number 


pulse 


phase If 


pulse 


phase ip 


1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 


Qf 

Qk+l 

Qf 

Qfc+i 

Qf-, 

Qf-i 

QIU 

QT 

Qf 

Qf 

Qfi 

Qf 

Qf 


-j-e 

77 - 89 + 1^ + a 
-f -37 + 30 
—a 

vr + 57 - 70 
TT + 27 - 20 + 6' + 2q; 


-87 + 100 

TT - 57 + 70 

2a 
-6 -a 
le + Aa 




Qf 

Qf-. 

Qf-. 

Qf^i 

Qf 

Qf 

Qk+l 

Qf 

Qk+l 

Qf-1 

Qf 

Qf 

Qf-. 

Qf 

Qf 


TT - 27 - 20 - §6' - a 
-f + 37-0 + ^ + a 

-f - 7 + 30 
-f -37 + 30 -0 + a 
-47 + 20 - 6* - 2a 
f-7-0-2^ + a 

¥ 




TT - 37 + 30 



\9 + a 






TABLE VI: Protocols for F-gates performing addition of the left addend qubit, k = L — 2. 

is mapped to another useful state according to the algorithm. If the transition is resonant 
the useful state is changed, if the transition is near- resonant, the useful state is not changed. 
The probability amplitudes of the useful states are not changed and are assumed to be real 
since the phases of the useful states are not important for the calculation of the probability 
errors: the unwanted states, generated from the useful states have random phases (see be- 
low), (b) Each existing unwanted state is mapped to another unwanted state according to 
the algorithm due to resonant transitions in the same manner as the useful states. (After 
generation, the unwanted states evolve as the useful states under the resonant action of 
the r/ pulses.) The existing phase of each unwanted state is not changed. Actually, this 
phase is not important because the phases of the amplitudes of the nonresonant transitions, 
which contribute to this unwanted state are assumed to be random, (c) Each useful state is 
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mapped to a finite number (which is equal or less than 2L) of unwanted states (nonresonant 
transitions). The details of this step are described in the rest of this Section, (d) The total 
probability of unwanted states is calculated and the probability amplitudes of the useful 
states are multiplied by a factor, close to unity, to satisfy the normalization condition. The 
total probability error is calculated as a sum of the probabilities of all unwanted states. 



A. Probability amplitudes for the nonresonant transitions 

The probability amplitudes for the nonresonant transitions were calculated in Ref . . If 
the state \i) is initially populated, Cj(to) = 1, then after application of one rectangular pulse 
with duration r and frequency z/, resonant or near-resonant with the transition frequency 
of the A;th spin, the probability amplitudes for the nonresonant transitions associated with 
flip of k'th spin {k' 7^ k) are 
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(18) 



Here the state \m) is related to the initial state \i) by flip of k'th spin; the state \n) is related 
to the initial state by flips of kth and k'th spin; o" = 1 if the k'th spin of the state \i) is in 
the state Uk' = and cr = — 1 if the A;'th spin of the state \i) is in the state n^' = 1. The 
other parameters are: 
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if rifc = then < 



A = E„ 



Er, 



D = Em - Ei - au 
+(A-(5)/2 



if rik = 1 then < 



?7 = —1 

S = E,- 

A = Ejn 

D = Er,- 



E,-u 
- E„ 



E, 



- V 
av 



+(A-(5)/2 



(19) 
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Here S is the detuning for the resonance or near-resonance transition \i) \j) associated 
with flip of the kth spin; A is the detuning for the resonance or near-resonance transition 
between the unwanted states \m) —>■ \n); is the state of the kth spin of the initial state 
The transition amplitudes for the non-resonant transitions in Eq. p8|l are characterized by 



the large detuning \D\ ^ \k — k'\6ui |l|,l2,y,l2|, which is approximately equal to the distance 
between the kth and k'th spins measured in the frequency units. 



B. Random phases for the unwanted states 

The analytical formula (jlSp and for the transition amplitudes of the nonresonant tran- 
sitions are derived for one pulse only. We now discuss how to use these equations for 
simulations of quantum protocols which consist of many rectangular pulses with different 
parameters. First assume that initially only one state is populated in the register of our 
computer. The first pulse of the protocol creates approximately 2L unwanted states as a 
result of the nonresonant transitions. The amplitudes of all nonresonant transitions, Cm, 
can be calculated using Eq. (fTHj) without integration of the Schrodinger equation. The 
probabilities of populations of these states are given by the square moduli Pm = |C.mP and 
the phases are equal to the arguments of Cm, <t>m = argC^. Since in Eq. (fT8|l to > t n/fl 
and D/Q 10"^ ^ 1 the values of (pm oscillate rapidly as functions of to and can be assumed 
to be random. In our simulations we assume random values for (pm in spite of the fact that 
the exact values are known, [(pm are equal to the arguments of the exponents containing 
Dtf) in Eq. (fTH|) .] We do this in order to make our method applicable to other possible 
quantum computer models where the probabilities of the nonresonant transitions can be 
estimated analytically, but the phases can not. On the other hand, as was shown above, the 
information about the exact values of these phases is not important. 



C. Linear accumulation of probability errors 

When the subsequent pulses of the protocol contribute to the same unwanted state, 
the corresponding probability amplitudes should be added to each other. The addition of 
complex numbers with the random phases can increase or decrease the modulus of the sum, 
depending on the difference between the phases of the addend numbers. Consequently, it is 
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possible, in principle, for the r/ pulse to decrease the total error. The result of the action 
of many pulses on one unwanted state can be presented as a multiple algebraic addition of 
complex numbers with random phases and the change of the probability amplitude of each 
unwanted state can be described by the random walk model in two-dimensional (complex) 
plane. Due to this model, the probability amplitudes of the unwanted states g row as K^rri ^ 
\fN , where is the number of pulses, and the total probability error is proportional to A^. 



VI. NUMERICAL RESULTS 



First, we simulate the full adder protocol for a small number of qubits using an exact 
numerical solution in order to calculate the phase errors and to test the quantum mar 
approach. For the relation between the numerical and physical parameters see Refs. 



m map 



A. Phase error 



The phase error, Pph, is defined as 

Pph = maXil^^G, - 501, (20) 

where 

Hg. = <PgXT)-'PbM. 0B.(O) = arctan |^g^'|°| , 0c,(T) = arctan |^g^'|^| , (21) 
where Im and Re stands, respectively, for the imaginary and real parts. The common phase 

1 M 

= JiT.i'pGAT) - <PB.m (22) 

i=l 

is close to the value —3a given by the protocol. The coefficients Cfi^(O), CdiT), and M 
are introduced in Eq. (fTUI) . The phase error in the FA(A) gate caused by nonresonant 
transitions is shown in Fig. 1 for different values of the addend number A for / = 5 addend 
qubits. (The total number of qubits is L = 11.) Initially, the superposition contained M = 4 
numbers: 7, 12, 16, and 27 with randomly chosen complex coefficients C_Bj(0) satisfying the 
normalization condition Y^fLi |Cb.(0)P = 1. The phase error in Fig. 1 is of the order or less 
than 1% of TT. This error can be decreased by increasing 5uj or decreasing VL Q]. The phase 
error increases as the number of qubits increases Q|. 
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Ql 

10 ^ 20 30 

FIG. 1: Phase error after implementation of the fuh adder for different values of the addend number 
A. K = 8[nf^ J/8, see Eq. ©], 5uj/n = 10^ {6uj w 1252J), / = 5 (L = 11). 

B. Test of the quantum map approach 

Since we used random phases for probability amplitudes generated by the unwanted 
transitions instead of their actual phases, we tested our approach by comparison with the 
exact numerical solution. In Fig. 2 we compare the probability errors calculated using 
quantum maps with the probability errors computed using the exact solution for L = 9 
qubits. The initial superposition contained the numbers 2, 5, 11, 12 with randomly chosen 
normalized coefficients (7^.(0) and A=6. The probability error for the exact case is defined 
as 

M 

P = T.\\CGXT)\'-\CBX0n (23) 
j=i 

while the probability error using quantum maps is defined as the sum of the probabilities of 
all unwanted states. As follows from Fig. 2, there is good correspondence between the data 
obtained using the exact numerical solution and the results obtained using the quantum 
maps. We observed a similar correspondence for the other initial conditions and other 
addend numbers A. 
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pulse number 



FIG. 2: The probability error as a function of the number of Q-pulses during implementation of the 
full adder obtained using exact solution and using quantum maps. K = 100 (J7 0.01), 5uj = 100, 
Z = 4 (L = 9). The data obtained using the quantum maps are averaged over 100 realization with 
different sets of the random phases. 




Pulse number 



FIG. 3: The probability error as a function of the number of Q-pulses during implementation of 
the full adder obtained using quantum maps. K = 100 {CI ^ 0.01), Su = 100. Number of addend 
qubits is I = 1000 and the total number of qubits is L = 2001. 
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Pulse number 



FIG. 4: The growth of the number of unwanted states generated by one useful state as a function 
of the number of Q-pulses implementing the full adder protocol. The parameters are the same as 
in Fig. 3. The data are averaged over 100 quantum map realizations. 

C. Modeling the full adder with 1000 addend qubits 

The growth of error with increasing the number of pulses is shown in Fig. 3. This figure is 
obtained using quantum maps with 20 randomly chosen initial numbers, Bi, with randomly 
chosen normalized coefficients, Cs.(O). The data are averaged over 20 different realizations 
with different sets of the random phases. The error bars are of the order of or less than 
the fine width. We also modeled the full adder for superposition of 1 and 100 states and 
obtained the same curve as in Fig. 3, so that the probability error appears to be independent 
of the number of states in the superposition. During modeling only the unwanted states with 
probabilities $, > 0.001(n/(5a;)^/M were taken into consideration. All other unwanted states 
with the smaller probabilities were neglected. Decreasing the value of C, did not significantly 
affect the value of the probability error. As discussed above, the probability error in Fig. 3 
grows linearly and is approximately equal to N{fl/Su!)'^, where is the number of pulses. 
In Fig. 4 we plot the number of unwanted states generated by one useful state as a function 
of the number of Q-pulses. If the initial superposition contains M useful states the number 
of unwanted states must be multiplied by M. 
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VII. CONCLUSION 



We have successfully demonstrated in this paper that a very efficient mapping procedure 
can closely approximate exact quantum dynamics of quantum computer. Our quantum 
map approach can also be generalized for use in other models, including models with time- 
dependent Hamiltonians. (Our Hamiltonian is time-independent in the rotating frame.) 
For a model with nearest neighbor constant interaction between the qubits, the result of 
the action of a pulse depends on the frequency difference between kth and k'th spin and 
on orientations of the following spins: kth spin with resonant or near-resonant frequency; 
{k — l)th and {k + l)th spins; k'ih spin with nonresonant frequency, whose flip creates 
an error; {k' — l)th and {k' + l)th spins. The number of these possible configurations is 
relatively small even for a computer with large number of qubits. Even in the case when 
the analytic solution is unknown, the results of the action of a single pulse on the finite 
number of possible spin configurations enumerated above can be calculated numerically to 
define the moduli for the transition amplitudes of one step of the discrete map. As shown 
in this paper the phases of the generated unwanted states can be assumed to be random. 
After identifying the maps, one can use them to simulate a whole quantum protocol. 

The quantum map method is similar to dynamical simulations of quantum systems with 



m 
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18l |. In the latter approach 



time-periodic Hamiltonians using the Floquet theorem 
one calculates the dynamics for one period of external field, T. Then one builds the evolution 
operator for one period, and then the evolution operator for one period can be used to find 
the state of the system at time mT, m = 2,3. . ., without integration of the Schrodinger 
equation for this time. The Floquet theorem is especially useful for finding the asymptotic 
behavior of the system, for m oo. Similarly, the quantum map approach is most useful 
when a quantum protocol contains a large number of pulses and when the number of qubits 
in the register of a quantum computer is large. 
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